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ABSTRACT 

When they first appear in the HR diagram, young stars rotate at a mere 10% of their break-up 
velocity. They must have lost most of the angular momentum initially contained in the parental 
cloud, the so-called angular momentum problem. We investigate here a new mechanism by 
which large amounts of angular momentum might be shed from young stellar systems, thus 
yielding slowly rotating young stars. Assuming that planets promptly form in circumstellar 
disks and rapidly migrate close to the central star, we investigate how the tidal and magnetic 
interactions between the protostar, its close-in planet(s), and the inner circumstellar disk can 
efficiently remove angular momentum from the central object. We find that neither the tidal 
torque nor the variety of magnetic torques acting between the star and the embedded planet are 
able to counteract the spin up torques due to accretion and contraction. Indeed, the former are 
orders of magnitude weaker than the latter beyond the corotation radius and are thus unable 
to prevent the young star from spinning up. We conclude that star-planet interaction in the 
early phases of stellar evolution does not appear as a viable alternative to magnetic star-disk 
coupling to understand the origin of the low angular momentum content of young stars. 

Key words: Stars; formation - Stars; rotation - Accretion, accretion disks - giant planet 
formation - Magnetohydrodynamics (MHD). 


1 INTRODUCTION 

It has long been known that the Sun and its siblings are character¬ 
ized by quite modest rotational velocities, of order of a few km/s at 
most (e.g. |Kraft|1970^ . This stands in sharp contrast to more mas¬ 
sive stars that exhibit much higher spin rates all over their evolu¬ 
tion. [Schatzm^ ([196^ first theorized that magnetized winds from 
solar-type stars would carry away a significant amount of angular 
momentum, thus braking the stars quite efficiently. Hence, regard¬ 
less of their initial velocity, stars with outer convective envelopes 
would end up as slow rotators on the main sequence. This expec¬ 
tation was later confirmed by |Skumanich| (1972) who showed that 
the rotational velocity of solar-type stars appeared to steadily de¬ 
crease on the main sequence, following the well-known Q. oc 
relationship. 

Extrapolating this relationship back in time to the early pre- 
main sequence (PMS), young stars were thus expected to be fast ro¬ 
tators at birth. More generally, as gravity dominates the late stages 
of protostellar gravitational collapse, newly-born stars were com¬ 
monly thought to rotate close to break-up velocity. Surprinsingly, 
the first measurements of spin rates for low-mass PMS stars, the 
so-called T Tauri stars (TTS), did not meet these expectations. On 
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the contrary, young stars were found to have only moderate rota¬ 
tional velocities, on average about 10 times that of the Sun’s, i.e., 
a mere 10% of their break-up velocity ( [Vogel & Kuhi|19^|Hart-| 
jmann et al.|1986[ [Bouvier et al.|1986| l. More than 30 years later, 

this aspect of the so-called ’’initial angular momentum problem” 
remains very much vivid. 


Several physical processes have been proposed to account for 
the slow rotation rates of young solar-type stars. They all rely on 
the magnetic interaction between the young star and its circustel- 
lar disk. At least 3 classes of models can be identified: X-winds 
d Slm et al.|1994^ , accretion-powered stellar winds jMatt & Pudrit^ 
[2005 I, and magnetospheric ejections ( jZanni & FetTeira|2013) . Ah 
these models investigate how the angular momentum flux between 
the star and its surrounding is modified by the magnetospheric in¬ 
teraction with the inner accretion disk. A brief outline of these mod¬ 
els and a discussion of their specific issues can be found injBouvierj 
jet al.j ( |2014| l. While a combination of processes may eventually 
provide strong enough braking torques on young stars to account 
for their low spin rates, none have definitely proved to be efficient 
enough. 


We investigate here an alternative process to account for the 
low angular momentum content of young stars. Specifically, we 
explore the flux of angular momentum being exchanged within a 
system consisting of a magnetic young stars surrounded by a close- 
in orbiting planet embedded in the inner circumstellar disk. We re- 
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view tidal and magnetic interactions between the protostar and a 
close-in planet to estimate whether spin angular momentum can 
be transfered from the central star to the planet’s orbital momen¬ 
tum and eventually from there to the disk hy gravitational interac¬ 
tion, thus effectively spinning down the central star. Previous papers 
have extensively explored star-disk, star-planet, and planet-disk in¬ 
teractions, usually with the aim to investigate the orbital evolution 
of close-in planets (e.g. |Zhang & Penev||2014| |Laine et al.|2008| 
ILaine & LinpOlIlpiMg et al.|2010||2012||Lai|2012| ). We build 

here from these earlier studies to explore the star-planet-disk inter¬ 
action in the specific framework of the initial angular momentum 
problem. Indeed, while addressing the issue of halting planet mi¬ 
gration near young stars through magnetic torques, |Fleck| j2008^ 
mentions that such interactions may be at least partly responsible 
for the slow rotation rates of PMS stars. We will show here that, for 
realistic sets of stellar and planetary parameters, magnetic as well 
as tidal torques seem actually unahle to prevent pre-main sequence 
stellar spin-up. 

In Section 2 we describe the general idea developed here, sum¬ 
marize the parameters of the system, and express the requirement 
for an effective protostellar spin down. In Section 3 we summarize 
the accelerating torques acting on the central star, namely accretion 
and contraction, that tend to spin it up. In Section 4, we explore 
both tidal torques and magnetic torques acting between the young 
star and the inner planet and review their efficiency in removing 
angular momentum from the central object. We discuss the quanti¬ 
tative torque estimates and the limits of the model in Section 5, and 
highlight our conclusions in Section 6. 



Figure 1. A sketch of the star-planet-inner disk interaction, with the main 
torques indicated: accretion, contraction, star to planet, planet to disk. 


Table 1. System parameters 


Star: M*= 1 Mq; i?*= 2 Sq; n*= 5 Hq; B* = \kG 

Disk/wind: Macc= IQ-^Moyr-*; M“’= 

Planet: Mpi= 10“^ Mq; Rpi=0.2Rq\ t„„j= 1 Myr 


2 STAR - PLANET - INNER DISK INTERACTION 
(SPIDI): A GENERAL ERAMEWORK 


The general idea that we develop here is to investigate whether a 
fraction of the spin angular momentum of the star can be transferred 
to the orbital angular momentum of a close-in planet, which in turn 
could lose its excess angular momentum to the disk through tidal 
interaction. In this way, the planet would act as a ’’lift” carrying an¬ 
gular momentum from the central star to the outer disk (cf. Fig.[TJ. 
We aim here at computing an equilibrium solution that would fulfill 
this requirement, and hence prevent the star from spinning up as it 
accretes and contracts. Alternatively, a more dynamical view of this 
process can be thought of, where the planet migrates successively 
inwards and outwards several times just outside of the corotation 
radius, thus transferring over time a significant fraction of the star’s 
angular momentum to the outer disk. 

In the following subsections we compute the torques acting 
in a system consisting of a young star, its circumstellar disk, and 
a close-in orbiting planet. The parameters of the system are listed 
in Table [T] To set orders of magnitude, we compute the specific 
angular momentum of the star, of the planet, jpi. 

The specific spin angular momentum of the central star is 


7 * 


— = 
M* 


( 1 ) 


with k^=0.205 for a completely convective star, which yields y* = 
5.8 10*^ cm^s“*. The specific orbital angular momentum of a non¬ 
eccentric planet orbiting at a distance a from the star is: 

jpi = (2) 

which is of order of a few 10** cm^s“* for an orbital semi-major 


axis of a few stellar radii. The specific angular momentum of the 
planetary orhit is 2 orders of magnitude larger than the specific an¬ 
gular momentum contained in the stellar spin. However, a Jupiter- 
mass planet orbiting at the corotation radius has nearly the same 
total angular momentum as the star, 

Jp, _ Mp,(GM,ay/^ _ Mp, 

for a = Tco = the break-up velocity fljr = 

, and the system parameters defined in Table[^ 

3 SPIN-UP TORQUES: YOUNG STARS OUGHT TO 
ROTATE EAST 

In this section, we describe torques that act to spin up young stars 
as they start their evolution on Hayashi tracks, namely the accretion 
and the contraction torques. 


3.1 The accretion torque 

A magnetic young star accreting from its circumstellar disk as¬ 
sumed to be in Keplerian rotation undergoes a spin-up torque ex¬ 
pressed as 

Tacc = MaccGQ.Kep{r,) = MacAGM , (4) 

where Macc is the mass accretion rate onto the star, r, is the inner 
disk truncation radius, Clfcep the Keplerian velocity in the disk, and 
M^, the stellar mass. 

For the system parameters listed in Tahle[T] the stellar magne¬ 
tosphere truncates the disk at a distance of 

r, ^ M-JPr'-^''' = 3.6 10" cm = SR^ = 2.51?*, 


(5) 
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where 6 * is the star’s magnetic field and the stellar radius 
( [Bessolaz et al.|2008| ). 

Alternatively, assuming that the truncation radius is located 
close to the corotation radius in the disk 

r, r,, = = 12R^ = bR^, ( 6 ) 


where is the star’s angular velocity. 

For these values of the truncation radius, the accretion torque 
amounts to T^cc^ 4.3-6.7 10^^ g cm^ s“^. The spin up timescale of 
the central star due to accretion is 




(7) 


with the stellar angular momentum 7*= 1.2 10^® g cm^ s“*. Using 
the system parameters listed in Table [T] the spin up timescale thus 
amounts to 7 Myr and reduces to 0.7 Myr for a protostellar mass 


accretion rate of 10 ^M^yr '. Flence, as shown by Flartmann & 


|Staufrer| ( |1989[ l, accretion from the circumstellar disk is expected 
to spin up the star to a significant fraction of its break-up velocity 
within a few Myr. 


3.2 The contraction torque 


As newly-born stars descend their Hayashi track, their radius de¬ 
creases and they eventually develop a radiative core. Both effects 
yield a reduction of the stellar moment of inertia and, if angular 
momentum is conserved, the star spins up. Taking into account only 
radius contraction during the first few Myr, the fully convective star 
can be described as a n = 3/2 polytrope. 

The potential energy of a poly trope of index n, is given by 


3 GMl 
5-n R^, 


( 8 ) 


and the luminosity of a star undergoing homologous contraction is 


T* 


(37 - 4) dEp„ 
3(7 - 1) dt 


= AnRio-r 


eff 


(9) 


The torque that should be applied to prevent the star from spin¬ 
ning up as it contracts down on the pre-main sequence on a Kelvin- 
Flelmholtz timescale is 


dl , dRt, 

= n*- = ( 10 ) 

dr dr 

and combining equations m with n = 3/2 and 7 = 5/3, one 
finally gets 


Tkh 



r 

GM, * 


( 11 ) 


With the reference parameters listed in Table[T]and assuming T^ff = 
4000 K, one finds Tfcn- 5 10^^ g cm^ s“^. Flence, stellar contraction 
is equivalent to an accelerating torque whose magnitude is similar 
to that of the accretion torque. Both contraction and accretion will 
thus equally act to spin the young star up. 


4 SPIN DOWN TORQUES: CAN THEY PREVENT FAST 
ROTATION? 

In this section, we investigate the torque a planet orbiting outside 
the corotation radius would exert onto the central star. Such a planet 
would extract angular momentum from star and migrate outwards. 
At the same time, the planet feels the tidal torque from the disk 
through Lindblad resonances and tends to migrate inwards. An 


equilibrium configuration could result if the torques acting on the 
planet balance outside the corotation radius, thus ensuring a con¬ 
tinuous transfer of angular momentum from the star to the planet 
and from the planet to the disk (cf. Fig.[TJ. If the outwards net flux 
of angular momentum can counteract the accretion and contraction 
torques, the star will thus be prevented from spinning up. Whether 
such an equilibrium configuration can be reached, and whether it 
is able to effectively spin the star down, depends on the nature and 
strength of the star-planet-disk interaction. We first discuss tidal ef¬ 
fects, and then turn to magnetic interactions. 


4.1 Tidal torques 

We investigate whether tidal torques are strong enough to allow the 
orbiting planet to extract spin angular momentum from the star and 
to pass it on to the outer disk. In order to spin the star down, the 
planet must be located beyond the corotation radius. The inward 
migration torque the disk-embedded planet experience must there¬ 
fore be balanced by the star-planet tidal torque beyond the corota¬ 
tion radius for the lift to operate. 

The migration torque exerted by the disk on the planet can be 
expressed as 

d7 rjl 1 I /T 1 /T dct j jj! 

Td.sk = -M,AGM^y/^a-'^^- -(12) 

df 2 dt T„ig 

where t„„j is the planet’s migration timescale in the inner disk. 

We may first verify that the magnitude of the migration torque 
outside the corotation radius would be able to compensate for the 
accretion and contraction torques. The condition can be written as 

\Tdisk(ti > r^oil > \Tacc + Tk{j\ — 2 \Tacc\ ■ (13) 


Using the above expressions for the accretion and migration 
torques, and assuming the inner disk is located at the corotation 
radius (cf. Eq.[^, this translates to 



2 t„ 


Mac 


(14) 


which, adopting parameters from Tableyields a > Area- Hence, 
if a giant planet can be hold off beyond the corotation radius, it 
may effectively extract enough angular momentum so as to com¬ 
pensate for the accretion and contraction torques, thus preventing 
the central star from spinning up. 

Is the tidal effect between the star and the planet strong enough 
to prevent the planet from migrating inwards of the corotation ra¬ 
dius? The tidal torque is given b >0 



where the migration rate of the planet due to the tidal interaction 
with the star is given by: 



= sign(fl* - Qpi) 


9 

2 



Mjy 

Q* 


(16) 


where 12 ,,/ is the planet orbital rotation rate, 2 * is the tidal dissipa¬ 
tion factor (e.g. |Zhang & Penev|2014^. The sign of the migration 


* The factor of 2 difference between Eq. |l5| above and Eq.(18) of |Chang] 
|et aL| j2010^ arises from their use of the apsidal motion constant which 
amounts to half the tidal Love number used to estimate the tidal quality 
factor Q,t(cf.[Mardling & Lin|200^. 
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depends on whether the planet orbits within or beyond the corota¬ 
tion radius. The migration induced by the disk is always inwards 


I da \ a 

I /disk Tra's 


(17) 


In order to slow the star down, the planet must orbit beyond the 
corotation radius. Hence, the planet’s inward migration must be 
stopped before it reaches the corotation radius. We must then have 


( 





dfl\ 

/disk 


(18) 


for r ^ all the way down to the corotation radius. Combining 
equations [T^ and [m we compute the critical semi-major axis at 
which the disk and stellar torques have the same amplitude on the 
planet (see also|Lin et al.|1996ll. 


/^crit ~ 


\ 22 * j 



(19) 


At a semi-major axis larger than the planet will migrate in¬ 
wards, while for distance smaller than the migration will be 
outwards. In order to ensure outwards migration at or beyond the 
corotation radius, as required to brake the central star, must 
thus to be larger than r„. Using parameters listed in Table [T] and 
2*= 10®, one finds acri,= 2.4 10** cm = O.Sr^. The mag¬ 

nitude of the tidal torque thus appears unable to prevent the planet 
from migrating inward the corotation radius, at which point it will 
contribute to spin up the star. Indeed, the tidal torque decreases very 
steeply with distance as seen from Eq.[T^ Its amplitude at the coro¬ 
tation radius is about 4 orders of magnitude weaker than required 
to balance the migration and accretion torques (cf. Fig. |^, a re¬ 
sult which is barely sensitive to the parameters of the system (cf. 
Eq.[T3. 


4.2 Magnetic torques 

4.2.1 Models and kind of magnetic interactions 

In this section, we aim at estimating orders of magnitude of the 
magnetic torques exerted by the planet on the star. To do so, we 
model the planet and the star as objects with a given radial profile of 
electrical conductivities, possibly generating magnetic fields. Nat¬ 
urally, this also models magnetic interactions between a moon and 
its planet, such as lo and Jupiter. The following is thus valid for 
star-planet or moon-planet magnetic interactions. 

In the simplest case considered by some authors (|Laine et al.| 
|2008[ |Chang et al.||2010l |2012^ , vacuum is assumed between the 
star and the planet. If the star and the planet both generate a dipolar 
field, the magnetic torque is simply given by the cross product of 
the magnetic moment with the magnetic field produced by the other 
dipole (which tends to align the magnetic moments). Torque can 
also come from dissipative effects. Typically, any time-variation of 
the magnetic field in an electrically conductive domain will gen¬ 
erate eddy currents, leading to Joule dissipation and torques. This 
is the so-called Transverse Electric (TE) mode ( |Laine et al.|2008l l, 
and the associated torques are e.g. studied by |Chang et al.| ( [2012[ l in 
the case of a planet orbiting in a tilted stellar magnetic field, which 
generates eddy currents in the conducting planet. 

Actually, the space between the star and the planet is rather 
filled by a stellar wind originating from the star, which can be 
roughly modeled as an electrically conducting fluid in motion. This 
has, in particular, two consequences: (i) the stellar magnetic field 


is adverted by the stellar wind, which is the so-called interplan¬ 
etary magnetic field (IMF), and decays thus less rapidly than in 
vacuum, and (ii) waves and currents can exist between the star and 
the planet allowing new kinds of magnetic interactions. The former 
point leads to a stronger TE mode in the planet, due to the higher 
time-varying magnetic field adverted by the stellar wind. The latter 
point requires to investigate the interaction between the planet and 
the surrounding magnetized flow of the stellar wind. 

The star-planet interaction via the stellar wind can be of var¬ 
ious kinds, depending of the planet’s velocity Vorb relative to the 
stellar wind (SW) local one v**, the speed Va,/ of the fastest wave 
in the stellar wind (i.e. the so-called fast magnetosonic wave), and 
the speed va of the shear (or intermediate) Alfven wave. First, if the 
(fast) Alfven Mach number Maj = llVorb “ VswII/va,/ is larger than 
1, the obstacle constituted by the planet generates a shock wave. In 
this so-called super-Alfvenic case, the flow is controlled upstream 
and disturbances are transmitted downstream. This is the case for 
Venus for instance, where the shock wave takes the usual form of 
a (detached) bow shock due to its blutf (i.e. non sharp-nosed) ge¬ 
ometry. In the case where the planet generates its magnetic field, 
like the Earth, this changes the apparent radius of the obstacle 
(aka the magnetosphere radius) for the stellar wind, leading to a 
enhanced coupling. This case of magnetic interaction between the 
stellar wind and the planet is called magnetospheric interaction by 
|Zarka| ( (2007l l. 


Second, we consider the sub-Alfvenic case (Maj < 1), fo¬ 
cusing first on the simple case where the planet is weakly mag¬ 
netized or unmagnetized. In this case, the planet motion generates 
shear (or intermediate) Alfven waves which propagate in the stel¬ 
lar wind and transport some energy from the planet to the star. 
The wave can then be reflected back, from the star to the planet, 
which leads to a planet-star interaction. Since the disturbance is 
only partially reflected, each subsequent reflection is of lower am¬ 
plitude, and results in a lesser change to the current system jCrary| 
|& BagenaT||1997^ . After several round-trip travel times, the sys¬ 
tem reaches a steady state, and a current loop is settled. In this 
case, the star-planet interaction can be modeled as a DC circuit 
( [Goldreich & Lynden-Bell|1969|l, which is c alled the unipolar in¬ 
ductor model or the TM mode ( Laine|20T3| l. In the other limiting 
case, the planet has moved away when the wave comes back and 
the planet’s interaction is thus decoupled from the star. This is the 
so-called Alfven wings (or Alfvenic interactions) model (e.g.|Saur| 
|et al.|20()4l >. These two models, sometimes presented as two differ¬ 
ent kinds of interaction, originates actually from the same phenom¬ 
ena and have been primarily applied to the lo-Jupiter interaction, 
but also to planets around magnetic dwarfs, ultra compact WD bi¬ 
naries, exoplanetary systems, etc. 

In the regime Maj < 1, we finally have the case where 
both the planet and the star generate a magnetic field. In the stel¬ 
lar system, the only example of such a case is the Ganymede- 
Jupiter interaction, which constitutes a textbook example of the 
expected plasma environment around close-in extrastellar planets 
( |S^[M4l l. This kind of interaction is often termed the dipolar 
interaction (e.g. |ZaA|2^. 


4.2.2 A generic torque formula 

Focusing on dissipative torques, one can link the dissipated power 
Pj with the star-planet dissipative torque T by 

/’<; = r||vp,-v«||/fl. 


( 20 ) 
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where a is the star-planet distance, Vpi is the planet velocity (v^/ 
Slpia) and v® is the advective velocity of the magnetic field at the 
planet orbit (i.e. the SW velocity v® = v®”, or v® = Cl^a in absence 
of SW flow). Dimensional arguments allow to write any torque ex¬ 
erted on the planet by the star as 


T = eaA,ff (21) 

where e is a dimensionless coefficient, a is the star-planet distance, 
A^f f is the typical surface area on which the star-planet magnetic 
coupling is effective, and P™ is the typical energy density of the 
coupling on the area (which can be seen as the ram pressure of 
the stellar wind when a stellar wind is considered). Focusing only 
on the magnetic coupling, this formula reduces to 

T = £aA,ff^ =eanRl,^^, (22) 

Po Pa 

where B^/f is a typical magnetic field value exerted on the area 
Agff, fio is the magnetic permeability of the vacuum, and = 
■fAfffJn is the typical radius of A^f f. 

Equation ( |22| requires the knowledge of e, B^ff and A^ff. The 
expression of e depends on the magnetic interaction we consider, 
which is discussed in section [4.2.3| The calculation of magnetic 
field Beff depends on the medium between the planet and the star. 
If this is vacuum, then Beff is simply deduced from the Gauss coef¬ 
ficients of the stellar magnetic field. If a stellar wind exists, B^f f is 
simply given by the model chosen for the stellar wind. Finally, Agff 
depends on the planet. In absence of any external conductive layer, 
Rohs = Rph which is thus a lower bound for Rots - However, a planet 
has usually an ionosphere due to photoionization of the upper at¬ 
mosphere (even without any own magnetic field), which leads to 
an induced magnetosphere above the ionosphere. This also means 
that the time-varying stellar magnetic field can be shielded by the 
ionosphere such that induced eddy currents in the planetary interior 
is small. Note that, as argued by |Chang et al.| ( [2012| l, atmosphere 
circulation may maintain this so-called exo-ionosphere on the per¬ 
manent nightside of a tidally locked planet. In this case, the radius 
Ri,i,s to consider is the exo-ionospheric radius, which is typically 
Piono = 1.1 - lARpi in the case of lo jZarka|[?007| l. Finally, the 
planet can also have a own magnetic field, either stored in its rocks, 
or generated by dynamo. In this case, the planet develops its own 
magnetosphere. The magnetopause radius Rmp is then the appar¬ 
ent radius R„i,s to consider, which can be estimated by the balance 
between the total ram pressure p^am coming from the star and the 
planetary magnetic field Bpi\ 


Rmp — Rn 


2 Pa Pn 


(23) 


where m depends on the planetary magnetic field geometry (m = 3 
for a mainly dipolar magnetic field, m = 4 for a mainly quadrupolar 
magnetic, field, etc.), and where 

1 1 

Pram = -^P l|Vp| - V™||^ -l- T — • (24) 

2 2 yUo 

Finally, R„ts is thus simply given by 


Robs — Riano^ Rmp)‘ (25) 

Particular forms of the generic formula \22) has already been con¬ 
sidered in previous works. For instance, neglecting Vpi in equation 
l |20| allows to recover equation 9 of |Zarka| ( |2007^ , where B^ff = 
is the IMF component perpendicular to the SW flow in the planet’s 
frame. Similarly, equation (9) of|Chang et al.|j2010) is recovered 


with equation {22\ , assuming no stellar wind and a dipolar stellar 
field, such that B^ff - B-^{Rtld)^, and considering = f2Rpi. 

As stated by |Zarka| p007| >, this general expression is simply 
the fraction e of the magnetic energy flux convected on the ob¬ 
stacle (aka the planet), and is expected to provide a correct order 
of magnitude whatever the interaction regime (unipolar or dipolar, 
super- or sub- Alfvenic) as long as the obstacle conductivity is not 
vanishingly small. The physics of the interaction is now hidden in 
the coefficient e, which is the difficult estimate we have to obtain. 


4.2.3 Application of the formula 


In this work, we aim at studying if the star can be slowed down by 
planetary magnetic torque. Since the super-Alfvenic regime only 
allows disturbances to be transmitted downstream, this requires to 
consider sub-Alfvenic regime such that a torque can be exerted by 
the planet on the star. Given that the stellar wind accelerates when 
flowing away from the star, this imposes to consider orbital distance 
a smaller than the so-called Alfven radius given by (jBessolaz et al.| 
|2008l l 


rA 


BiRf 

-— = 12R. 

2GM,MlJ 


(26) 


distance at which the SW flow becomes super- Alfvenicj^ We have 
thus to estimate torques in four cases, the TE mode, the Alfven 
wings model), the dipolar interaction, and the unipolar inductor 
model. 


Torque due to the TE mode: In this case, the torque is associ¬ 
ated with the Joule dissipation in the interior of the planet, origi¬ 
nating from the stellar magnetic field oscillating at the rate oj. Ow¬ 
ing to the diffusive nature of the problem, one can thus estimate 
the Ohmic dissipation by ( |Cam p bell|1997| |1983[|Laine et al. |2008[ 
|Chang et al.|2012} : 

Rob,nic=^OJ^, (27) 

2pa 


where 'V is the volume where eddy currents take place. This is thus 
typically a thin shell of thickness 6, which gives 'V = A^ff S, i.e. 
'V ~ 47rBp,5 if it takes place at the surface of an electrically con¬ 
ductive planet. Noting rj the typical planetary electrical diffusivity, 
S - ( 277 / 01 ) is then the skin depth for our magnetic induction prob¬ 
lem. 

Einally, considering that w ~ Op/, an upper bound of the 
torque is obtained for'T' = 4;r/3 and we will thus use equation 

using Boi, = Bp, and 


Alfven wings: As shown by |Drell et al.| ( [T965^ , the radiation of 
two pure linear Alfven wings leads (when converted in S.I. units) 
to a torque given by formula j22| l, with the coefficient (see also |Lai| 

[M 2 l > 


e = 2MA, 


(29) 


^ Using the expression of the Alfven radius from |Matt & Pudrit^j2008) ’s 
stellar wind models (their Eq.l4) would result in twice as large an estimate. 
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Figure 2. Magnitude of the torques (in accretion torque units). Considering the left vertical axis, the curves show the magnetic torque of the TE mode (eq. |28| 
and Robs — ^pl), of tho Alfven wings (eq. |31| and Robs - the dipolar interaction (e = 1, eq.j^for a dipolar planetary surface magnetic field of lOG), the 
tidal torque, the magnetic torque associated with the unipolar inductor model (eq. |32| with a planetary conductivity of CTpi = 10“^ S/m.), and the migration 
torque (r^/g = 10^ yr), shown by a dotted line. From left to right, the solid lines are respectively the corotation and the Alfven radii, whereas the vertical dashed 
line shows r^, the upper radius of the favorable DC circuit closure regime. We have also plotted, for comparison, the empirical formula obtained by Strugarek 
et al. (2014) from 2.5D numerical simulations. Note also the change of slope noticeable at a distance of 5 in some curves, which is due to the stellar wind 
becoming unable to prevent the magnetic planet to develop an extended magnetosphere. Only the magnitude of the torques is plotted, not their sign, which is 
the reason why the tidal torque (cf. Eq. [ID does not change sign at the corotation radius in this plot. 


where M .4 = ||Vorb “ VswII/va is the shear (or intermediate) Alfven 
wave Mach number, with va = s/pjJo is the shear (or intermedi¬ 
ate) Alfven wave velocity and p the stellar wind density. The factor 
2 in equation l |29| > naturally comes from the existence of two Alfven 
wings. The work of|Drell et al.|jl965|l has then been extended by 
|Neubauer|jl980t to the fully non-linear sub-Alfvenic situation, also 
including flow v„rb - Ssw which is not perpendicular to the back¬ 
ground field B. Following |Zarka| p007| l, and taking into account 
the existence of 2 wing, this leads to the extended formula 


where the flow is perpendicular to the magnetic field (Vrei • B = 0) 
leads then to ( |Zarka|2007| l 


1 


(31) 


for the Alfven wings interaction. Note that, in this case, Rgi,s is of 
the order R^bs ~ Rph even in the presence of an exo-ionosphere. 
Note also that Alfven waves cannot be radiated when r > ta, and 
the Alfven wings interaction is thus only valid for r < va- 


2Ma 

e = - 

IIMAVrel -l-BII 


(30) 


with the two unit vectors Vrei and B, directing respectively the rel¬ 
ative velocity v„rb - Vsw and the (background) magnetic field B. 
Note that we naturally recover equation l |29| > for Ma « 1. Focusing 
on orders of magnitude, the factor 2 can be dropped, and the case 


Dipolar or magnetospheric interaction: As argued by |Zarka| 
( |2007) , the torques associated to these interactions can also be mod¬ 
eled by formula \22\ , using e = 0.1 - 1 (typically, e ~ 0.3 in the 
sub-Alfvenic Ganymede-Jupiter interaction). Note that e = 1 is in 
agreement with equation l |31| l, which gives e ~ 1 for Ma ^ 1. The 
important feature of this kind of interaction is the presence of a own 
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planetary magnetic field, which can lead to large magnetic interac¬ 
tions (due to large RchP. 


Unipolar inductor model (DC circuit or TM mode): In this case, 
first studied by |Goldreich & Lynden-Ben| ( |1969| l, the Alfven waves 
round trips between the star and the planet allow to reach a steady 
state, which can be modeled as a DC circuit. In the planetary frame, 
the electrical field E = v x B generated by the planetary velocity v 
allows indeed to close a DC circuit through the surrounding elec¬ 
trically conductive medium. An electrical current is thus flowing 
through the electrical resistances of the planet, Rp, of the stellar 
footprint of the planet, '/?*, and of the two flux tubes crossing the 
stellar wind, IRf. 

The difficult estimation of the three electrical resistances is 
performed here by following |Laine & Lin| j2012| >; |Laine| ( |2013| l, 
who neglect ^ ??p, ??*. According to the formula obtained by 
|Laine & Lin| ( |2012t ; |Laine|j2013^ , we thus end up with 

8^ |Dp/— D*| 


One can distinguish different regimes of closures of the DC circuit 
( |Laine|2013[ ). In the so-called favorable DC circuit closure regime, 
the Alfven wave have the time to perform the round-trip planet- 
star-planet, which actually imposes here that the planet is on an 
orbit below a certain distance In our model, r, ^ 8/?*. Analytical 
estimates of r,, are derived in Appendix A. 

Note that the DC circuit is naturally opened for r > r^. Note 
also that it exists an upper limit to the magnetic interaction torque 
generated using the unipolar inductor model ( |Lai||2012| ). Indeed, 
when the circuit resistance is too small, the large current flow 
severely twists the magnetic flux tube connecting the two binary 
components, leading to the breakdown of the circuit. According to 
|Lai|(|2012ll, this breakdown appears when 


4yUo a|n,,f - n*l _ £ 

n Kp+Kt 2 


(33) 


which limits 6 to e = 2 (corresponding to ^ = 1 in|Lai|20T2ll. 


4.2.4 Results 

In order to calculate the various torques estimates give above, a 
stellar wind model is required. In this work, we have considered 
two stellar wind models. One is adapted from |ZarS| ( |2007[ l, as de¬ 
tailed in Appendix B, and the other one is described in Lovelace 
et al. (2008). Both give similar results as far as magnetic torques 
are concerned (cf. Appendix B). Using the stellar wind model ad- 
pated from |Zarka| ( |2007| , the results are summarized in Figure]^ 
and compared with the tidal torque. Note that the magnetic torque 
associated with the DC circuit depends on the chosen planetary 
conductivity. It turns out that above CTpi = 10“^ S/m (as used in 
Fig.m, the resistance is very small > 1), and the torque is thus 
maximum, given by e = 2. 

Fig. clearly illustrates the main result of this study. None 
of the planet-induced torques acting on the central star, be it tidal 
or magnetic, is able to balance the accretion torque beyond the 
corotation radius. Hence, at least with the parameters adopted here, 
the star-planet interaction cannot prevent the spin up of the central 
star as it accretes from its circumstellar disk and contracts down 
its Hayashi track. Indeed, even the most powerful torques, corre¬ 
sponding to the dipolar interaction and the unipolar inductor model, 
fail by several orders of magnitude to match the magnitude of the 


migration torque beyond the corotation radius. We have also com¬ 
pared these results to |Strugarek et al.|p014| l 2.5D numerical sim¬ 
ulations and find that their parametric torque formulation predicts 
order of magnitude torques similar as the DC and dipolar cases 
investigated here (see Appendix|^for a short account of their para¬ 
metric torque formulation). 

One may wonder why the magnetic coupling seems to be effi¬ 
cient enough in the similar study of |Fleck| ( |2008^ , in contrast to our 
results. This is partly due to the fact that he prescribes the cross- 
sectional area A^ff of the magnetic flux tube linking a planet to its 
host star to be in the range a = \ — 10% of the stellar surface. Here, 
we determine A^ff self-consistently and show that it is of order of 
magnitude of where R„i,s ~ Rpi when the planet is close from 
the star, and thus a x kR^^iIAttR^ x 0.2% with our values. Also, 
|Fl^j2008t assumed a stronger magnetic field (2 kC) and a larger 
stellar radius (3/?0) than we adopted here, which results in a much 
stronger torque (cf. his Eq. 14). Indeed, a 3 kC field associated to a 
4 Rq protostar would be required to obtain a torque comparable to 
the accretion and migration torques. 


5 DISCUSSION 

We have considered the various torques acting in a system consist¬ 
ing of a contracting pre-main sequence star still accreting from its 
circumstellar disk into which a newly-formed Jupiter-mass planet 
is embedded on a short-period orbit. We have shown that a balance 
between accelerating torques acting onto the star, due to accretion 
and contraction, and decelerating ones, due to the star-planet inter¬ 
action, cannot be reached under the conditions we explored here. 
For all cases we investigated, the accretion torque exceeds any plan¬ 
etary torque acting on the central star by orders of magnitude. 

The main limitation of the scenario outlined here lies in the 
adopted wind model. The model is an extrapolation from a solar- 
type wind properties, which is not necessarily adapted to the wind 
topology of young stars. Firstly, the topology of the magnetic field 
of young stars might be quite different from that of a mature solar- 
type stars (e.g. |Donati & Landstreet|2009| l , although current mea¬ 
surements suggest that fully convective PMS stars at the start of 
their Hayashi tracks host strong dipolar fields, of order of a few kC 


(|Gregory et al.|2012| 

Johnstone et at. 2014b. Secondly, as shown 

by|Zanni & Ferreira| 

20091, the structure of the stellar magneto- 


sphere might be significantly impacted by its interaction with the 
circumstellar disk. In particular, an initially dipolar magnetic field 
may evolve into a much more complex and dynamical topology. 
Whether this would significantly modify the torques associated to 
the magnetic star-planet interaction is yet unclear. 

Another issue of course is whether the framework proposed 
here is applicable to the vast majority of young stars that rotate 
slowly as they appear in the HR diagram. This would require that 
planet formation be not only a common occurrence, but also that 
it is fast enough to impact the earliest stages of stellar evolution, 
and that it is quite dynamic so as to frequently send massive plan¬ 
ets on inner orbits. Whether all these conditions are met around 
protostars is yet unknown. The recent ALMA image of the numer¬ 
ous gaps in the disk of the protostar HL Tau suggests that multiple 
planet formation may indeed proceed quite rapidly after protostel¬ 
lar collapse. The significant fraction of hot Jupiters found to revolve 
around their host star on a non coplanar orbit further suggests that 
dynamical interaction between forming planets in the protostellar 
disk may be efficient to scatter massive planets close to their par¬ 
ent star jFabrycky & Tremaine|2007[|Triaud et al.|2010[ [Morton &| 
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|Johnson|201 l[[Albrecht et al.|2012[|Crida & Batygin|2014[ ). Even¬ 
tually, Hot Jupiters (HJs) are deemed to coalesce with the central 
star after a few hundred million years jLevrard e t al.| 2009||Bolmont| 
|et al.|20l^ . Hence, the scarcity of HJs around mature solar-type 
stars does not necessarily exclude that they are much more common 
during early stellar evolution (e.g.,|Teitler & Kdnigl|2014||Mulders| 
let al.|2015| >. 

Finally, we wish to emphasize that the preliminary results pre¬ 
sented here have to be extended in order to explore a much larger 
parameter space. Table[T]summarizes typical parameters for young 
accreting stars that we adopted here for the torque computations. 
However, classical T Tauri stars at a given mass and age exhibit a 
wide range of spin rates (|Bouvier et al.|2014|(, mass accretion rates 
(|Venuti et al.|2014) , magnetic field strength and topology ( |Donati| 
|et al.|2Q13^ . Henre, the torque estimates provided here may vary 
significantly from one system to the next. Moreover, during the em¬ 
bedded phase, other sets of parameters could be relevant, implying 
for instance a larger radius, and possibly stronger magnetic fields 


(e.g., 

Grosso et al. 

1997 

|Tsuboi et aL|2000| and faster rotation rates 

(e.g., 

Covey et al. 

2005 

L Whether the planetary lift scenario could 


be instrumental in the protostellar embedded phase to produce the 
low angular momentum content of revealed young stars remains to 
be ascertain. Also, more complex scenarios than that described here 
could be envisioned. For instance, the planetary migration torque 
could be significantly reduced, and the migration timescale consid¬ 
erably lengthened, by tidal resonances occurring between several 
inner planets embedded in the disk, or by considering a steeper sur¬ 
face density profile in the disk close to the magnetospheric cavity. 
In the latter case, the enhanced outward flux of angular momen¬ 
tum between the inner disk material and the close-in planet would 
result in sub-Keplerian disk rotation close to the truncation radius, 
thus potentially reducing the accretion torque onto the central star. 


6 CONCLUSION 

As an alternative, or a complement, to star-dik interaction mod¬ 
els that attempt to account for the low spin rates of newly-bom 
stars, we explored here the tidal and magnetic interactions between 
a young magnetic star and a proto-hot Jupiter embedded in the in¬ 
ner disk. We investigated whether such a close-in embedded planet 
could extract enough angular momentum from the central star to 
counteract both the accretion and contraction torques, thus leaving 
the star at constant angular velocity as it evolves on its Hayashi 
track. While a lot remains to be done in exploring the parame¬ 
ter space relevant to such systems, it appears that the decelerating 
torques acting on the central star are orders of magnitude too weak 
to counterbalance the spin up torques. Hence, even though plan¬ 
etary systems may promptly form around embedded young stars, 
as suggested by the recent ALMA image of HL Tau’s disl0 star- 
planet interaction may not eventually prove a viable alternative to 
magnetic star-disk interaction to understand the origin of the low 
angular momentum content of young stellar objects. 


APPENDIX A: ANALYTICAL ESTIMATES OE Ry EOR 
THE UNIPOLAR INDUCTOR MODEL 

Based on the work of |Laine| ( |2013t , we have been able to obtain 
analytical estimates of r^,, the transitional radius between favorable 


cf. http://www.eso.org/public/unitedkingdom/news/esol436/ 


and unfavorable DC circuit closure regimes (see. Fq. |32^ . Indeed, 
noting to = 2J?p;/||Vpi - v®”|| the advection time of the unperturbed 
plasma to cross the planet’s diameter, and tA,FT the Alfven wave 
travel time between the planet and the top of the footprint, the con¬ 
dition for favorable closure regime reads then to > Estimat¬ 

ing Ia.ft by the time Ia.fto it takes for the Alfven wave to travel 
along the flux tube until the stellar chromosphere, we obtain the 
following asymptotic estimates: 


ry 


RA 


4n E, 


H/4 


(Al) 


for slow rotators (Clpi « D* at r = r„), and 


= R^, 


9BlRl,R. 47r E, 


G/H* ^lo M'! 


1/5 


(A2) 


for fast rotators (Clpi » D* at r = r,,). 


APPENDIX B: STELLAR WIND MODEL ADAPTED 
EROM ZARKA (2007) 


The star we consider being not very different from the Sun, we 
have chosen to slightly adapt the solar wind model proposed by 
|Zarka| j2007^ . First, let remind roughly the various usual scaling 
laws at play in the solar wind, before detailing the exact model 
used in this work. At a certain distance r to the Sun (beyond a few 
Sun radius), the radial SW velocity can be considered as constant, 
and thus, mass and magnetic flux conservation give a SW density 
N and a radial field S,. decaying as 1/r^, whereas other conserved 
quantities (e.g. angular momentum) give an azimuthal velocity 
and magnetic field which decays as 1/r (e.g. jWeber & Davisj 
1 1 967[[Belcher & MacGregor| 1976^ . Close from the Sun, is rather 
dipolar and decays thus as 1 jr^ (which leads to a decay in 1/r^ of 
B^, see eq.|B2|(. 

The model proposed by |Zarka|p007| l for the Sun and the solar 
wind is (in spherical coordinates) 


Br 

B^ 

Bg 


{rlR.f y ^ (r/RA^ 

Q. r 

B, -, 

IIVpl - v®'*|| 

0 


(Bl) 

(B2) 

(B3) 


with the constants 4 < f < 9 (here taken equal to f = 6.5) and 
1.5x 10-^ T, and 


/V,(cm-3) = tti (fl/J?*)-‘5 -H a2 (alR,r'>/^ + 0:3 (a/R,)-^ (B4) 


withtti X! 3x 10* cm“*, Q '2 « 4x 10® cm“*, and 0:3 x: 2.3x 10® cm“®, 
which gives the SW density p = 1.1 nip Ne (with the proton mass 
nip X 1.673 X 10“^^ kg). In this model, the solar wind velocity Vsw 
can be described by 


S'!!' 

; 

CO 


1 - exp 



0 , 

0 , 


(B5) 

(B6) 

(B7) 


with V™' = 385 km/s. Note that the values given by this expression 
of V™ are naturally quite close from the ones derived from p by 
mass conservation (i.e. using pv™ = cst). 

We have adapted this model to our star by multiplying by 
M'^IR-l: the magnetic field by B*, and Vsw by a/M^JR^, 
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APPENDIX C: RESULTS OE STRUGAREK ET AL. (2014) 



Figure Bl. Comparison of the stellar wind speed V = ||Vsw “ Vorbll in the 
planetary frame for the model adapted from Zarka (2007), labeled Z07, and 
Lovelace et al. (2008) labeled LOS. The planetary orbital speed has 
also been plotted. 



Figure B2. Radial and azimutal magnetic field of the stellar wind for the 
model adapted from Zarka (2007), with the typical slopes (stellar magnetic 
field rescaled by its surface value of 1 kG). 


all these quantities being expressed in solar units (which allows 
to recover the initial solar wind model of |Zarka|2007^ , i.e. Rq = 
6.96 X 10* m, Mo = 1.99 x 10“ kg, In/Oe = 27 d, M"' = 1.6 x 
10“'^^ and Sq = 10 G (these two latter values are actually 

obtained as outputs from the solar wind model of |Zarka|2007r 

Note finally that theoretical models, such as the one of|W] 
|ber & Davis| ( |1967} , rather advocates an azimuthal magnetic field 
given by = 6r(v™ - as well as non-zero azimuthal 

velocity, typically given by co-rotation (vj"' = G*r) for r < ta and 
V™' = for r > ta- We have checked that this does not 


According to Strugarek et al. (2014) , the torque applied by the 
planet to the star is given by (see their Eq. 16) 

-5 


r = Kr^.c\^+b 


eo-®\ alR 


(Cl) 


where t„ = where c, b, p, t, 0 and i are given by their 

table 4, values obtained for a/7?* = 3, and where the last factor on 
the right hand side comes from their proposed scaling in (see 
their Eq. 28). We have added a supplementary factor K to remind 
that their values are all obtained for a planetary radius RpilRt = 
0.1, and the extrapolation to other planetary radius requires thus a 

(7?,,//7?*)-/O.P, since 


correction K. We propose for instance K 

~ ~ 


eq. 22 gives a torque in and eq. 23 gives Robs Rpi 
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